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We study the evolution of non-periodic cosmic string loops containing Y-junctions, such as may 
form during the evolution of a network of (p, q) cosmic superstrings. We set up and solve the 
Nambu-Goto equations of motion for a loop with junctions, focusing attention on a specific static 
and planar initial loop configuration. After a given time, the junctions collide and the Nambu-Goto 
description breaks down. We also study the same loop configuration in a U(l)xU(l) field theory 
model that allows composite vortices with corresponding Y-junctions. We show that the field theory 
and Nambu-Goto evolution are remarkably similar until the collision time. However, in the field 
theory evolution a new phenomenon occurs: the composite vortices can unzip, producing in the 
process new Y-junctions, whose separation may grow significantly, destabilizing the configuration. 
In particular, an initial loop with two Y-junctions may evolve to a configuration with six Y-junctions 
(all distant from each other). Setting up this new configuration as an initial condition for Nambu 
Goto strings, we solve for its evolution and establish conditions under which it is stable to the decay 
mode seen in the field theory case. Remarkably, the condition closely matches that seen in the field 
theory simulations, and is expressed in terms of simple parameters of the Nambu-Goto system. This 
implies that there is an easy way to understand the instability in terms of which region of parameter 
space leads to stable or unstable unzippings. 



I. INTRODUCTION 

Cosmic strings have long been known to arise naturally 
in a wide class of field theories [l], [3i S 01 j and recent work 
has indicated that they may also arise in superstring/M- 
theory l5|,j6[, for example through models of brane in- 
flation 0, la, Q- Back in the 1980's fundamental strings 
and cosmic strings were considered as separate. Obser- 
vational constraints on the allowed cosmic string tension 
had to be violated by fundamental strings in order for 
them to unify the gauge couplings at high energy p^ . 
Moreover, whereas cosmic strings could be of cosmologi- 
cal scales, fundamental type I and heterotic strings were 
known to be unstable and their natural size was of order 
the Planck scale [Toj . 

Things changed, however, with the second string rev- 
olution in which the importance of higher dimensional 
brane solutions was understood, including new types of 
strings, and it has since been realized that the extra spa- 
tial dimensions could be lar ge u sing warped compactifi- 
cations (for reviews see [TTI Il2|'). As a consequence, it 
turned out that the effective string tension of the fun- 
damental strings could be reduced, and depending on 
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the compactification scheme it became possible to have 
meta-stable fundamental strings of cosmological length 

This has opened up the fascinating possibility that this 
new class of cosmic strings, formed after a period of infla- 
tion, could provide a unique window on string theory, via 
measurements of their imprint on the cosmic microwave 
background radiati_on 16 , [H, [H, [13, HH , their lensing 

[or through their production 



of distant galaxies [22_ 
of gravitational waves and massive particles 

MM- 

When two cosmic strings intersect, the traditional lore 
dictates that they intercommute (swap partners), which 
in the case of a string looping back and intersecting it- 
self results in the loop being cut off. For traditional 
fleld theory strings (essentially Abelian Higgs strings 
close to the Bogolmonyi limit) this intercommutation oc- 
curs with effectively unit probability at each intersection 
jsol . Isil . [3^ . m, [11] . However in the cosmic superstring 
case the intercommutation probability would be reduced 
due to the presence of extra dimensions and a small string 
coupling constant [H, [H, [s^, HI] . 

Additionally, one of the properties of cosmic super- 
strings is that there can be two different types of string, 
along with stable composites of the two. That is, in addi- 
tion to fundamental (F) strings, there are Dirchelet (D) 
strings and also {p,q) composites of p F-strings and q 
D-strings. Where composites split up into more basic el- 
ements, there will be Y-shaped junctions, something not 
present in the cosmic string scenario that is convention- 
ally considered. This can lead to very different dynamics 
when strings collide and a crucial outstanding question 
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is how these would affect the properties of the string net- 
work and what the observational consequences would be? 

However, this question is not restricted to cosmic su- 
perstrings. Even the simplest model of gauge strings (the 
Abelian Higgs model) can provide stable composites if 
the gauge coupling is set sufficiently high. Indeed, for 
par t of this article we will use the double U(l) model 
[39| to explore the properties of the Y-junctions that it 
permits, thereby increasing the applicability of our in- 
vestigations, although it should be noted that the mass 
spectrum of the bound states is a specific prediction for 
cosmic superstrings that is not reproduced by these field 
theoretic models. 

The focus of this article is simple, a detailed compar- 
ison of the two approaches that have been adopted to 
describe the dynamics of strings containing Y-junctions. 
The first, developed by Copcland, Steer and Kibble in 
[13, El (CKS), is based on modifying the Nambu-Goto 
action to include junctions, and the second is to describe 
the strings as composite objects in terms of an underly- 
ing classical field theory that will allow for the formation 
of junctions [3§],[43l-[4§|. Such a comparison has a well 
established history for ordinary cosmic strings, where it 
has been shown that the Nambu-Goto action for a rela- 
tivistic ID string is an excellent approximation to the dy- 
namics of widely separated gauge cosmic strings as long 
as the curvature is on scales far greater than the mi- 
croscopic string width. Once that regime breaks down, 
the Nambu-Goto approximation is no longer such a good 
approximation to the full field theory (for a review see 
@j How does the comparison fare for these more 

complicated arrangements representing strings with junc- 
tions, and full loop configurations? Recently Bevis and 
Saffin showed that the results obtained by CKS for 
the collision of straight strings are consistent with those 
found using field theory simulations of the double U(l) 
model, once the composite region grows to be much larger 
than the string width. This is quite remarkable given the 
fact that the conditions for Nambu-Goto applicability are 
breeched in the region of the Y-junction at all times (a 
related result was also shown for a different field theory 
model in [sot). 

In the following, we employ the CKS approach in nu- 
merical simulations for the first time, exploring scenarios 
not available through analytical means and in doing so 
we provide a quite different comparison with field the- 
ory simulations to that previously performed. Of par- 
ticular note, we are able to address directly questions 
such as how closely the two approaches match, where 
the deviations arise and whether we can safely adopt the 
Nambu-Goto formalism to describe the evolution of loops 
containing Y-junctions. The answers to these questions, 
and in particular the last question is important. It is 
far easier to evolve a loop based on a modified Nambu- 
Goto action than on a full field theory action, hence it is 
possible to analyse a broader class of configurations and 
obtain more believable statistics concerning the distribu- 
tion of cusps and kinks on these loops - which in tm'n 




FIG. 1: An example of a loop configuration with multiple 
junctions: the butterfly configuration with a central string of 
tension /io and two arc strings of tension fj,i. The basic but- 
terfly conflguration has just two junctions, but we find that 
under certain situations these can decompose, as indicated 
by the magnified region, with a single junction splitting into 
three junctions that then continue to separate. 



is important if we want to establish gravitational wave 
templates for these objects. Motivated by an instability 
we see in our field theory simulations, we also study it in 
the context of the Nambu-Goto approach, both analyt- 
ically and numerically. This then allows us to establish 
the criteria in terms of a simple angular parameter, un- 
der which the junctions are stable to unzipping. This 
in turn will be important in the future when modeling a 
network of strings with junctions numerically using the 
Nambu-Goto equations. 

For simplicity, and to allow a direct comparison, we 
concentrate in this paper on a particular type of loop 
configuration, shown in Fig. [TJ It resembles a butterfiy 
and looks somewhat like the result of two co-planar loops 
colliding. We stress that this is not intended to represent 
a likely cosmological scenario, but presents an interest- 
ing case with which we may explore the properties of 
Y-junctions, including a new feature, their stability to 
decomposition into three separate Y-junctions. 

In sectionini we discuss the Nambu-Goto approach and 
describe our numerical technique, for loops with just two 
junctions and also loops with multiple junctions. This is 
followed in section Hill with a description of the field the- 
oretic simulations, and then in section HVl we present our 
results for each of the approaches including the identifica- 
tion of an instability in the field theory case that is then 
seen for the first time to also be present in the Nambu- 
Goto case, once suitable initial conditions are chosen. We 
conclude in section |V] and finish with three appendices 
giving detailed derivations of some of the results used in 
the main paper. 
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II. NAMBU-GOTO SIMULATIONS 
A. Equations of motion 

In this section we set up the Nambu-Goto equations 
of motion for a string loop with J junctions, generahsing 
the aforementioned CKS approach that studied the dy- 
namics of one junction of three strings. For example, a 
loop with a few junctions is shown in Fig. [1] where the 
colours differentiate the different charges and tensions of 
the strings. 



1. Case of two junctions 

For simplicity, we first derive the equations of motion 
for a loop with two junctions, which wc label by the in- 
dex J = {A,B), and three strings. We parameterize the 
position of the ith string (i = 1, 2, 3) as 



(1) 



where r and ct^ are the world-sheet coordinates (note that 
T is chosen to be the same for all three strings). The 
induced metric on the world-sheet for string i is 



Tab 



(2) 



where a — (t, cTj) and 77^,^ is the 4-dimcnsional Minkowski 
metric. Below, a dot/dash denotes a derivative with re- 
spect to r/ai respectively. The values of the world-sheet 
coordinate a at the junction are denoted by sf and are 
generally r-dependent. We are free to choose the direc- 
tion of increasing a on each string with respect to a given 
junction. In the case of two junctions, we may take a to 
increase from junction A to junction B so that 



\r) < a. < sfir). 



(3) 



For multiple junctions, as discussed below, this choice 
cannot always be made. The positions of the junctions 
are 



X^j{T)=xf{r,sUr)) 



for all i. 



(4) 



During the evolution of the loop, it may occur that 
for one of the strings and at a certain time Tc, s^(tc) = 
sf (tc). In that case, the physical length of the string 
between the junctions vanishes, and the junctions collide. 
More generally, junctions A and B collide when X^(tc) = 
Xg{Tc) which may occur even if sf{Tc) 7^ sf{Tc) and 
the length in string i does not vanish. Since in string- 
theory the outcome of the collision of junctions is not well 
understood, we will stop our NG simulations of loops at 
the collision time. The field-theory simulations presented 
in Sec. IIIII can, of course, go beyond this time. 

In the absence of background fluxes and after dilaton 
stabilisation, the dynamics of a single infinite {p, q)-stTmg 



in flat spacetime is given by the Dirac Born Infeld (DBI) 
action |5l| 



DBI 



Hab - 



(5) 



where p. = \q\/{gsX) is the tension of q coincident D- 
strings, A = 2'Ka' , with a' the Regge-slope parameter, 
and gs is the perturbative string coupling. Fab is the 
electromagnetic tensor on the string world-sheet, and the 
electric flux density is the momentum conjugate to the 
electric field p = dLoBi/dFra- The dynamics of three 
semi-infinite (p, g)-strings meeting at a junction was dis- 
cussed in (sTj where it was shown that the resulting equa- 
tions of motion arc exactly equivalent to those obtained 
by using the Nambu-Goto action for each string, provided 
the ith string tension in the Nambu-Goto action is taken 
to be given by 




(6) 

and one imposes charge conservation at the junction 

i i 

For this reason, and for simplicity, we assume that the 
dynamics of each individual segment of string is deter- 
mined by the Nambu-Goto action. 
In the conformal gauge 



(8) 



the Nambu-Goto action for the three strings of tensions 
/ii joined by two junctions is 

S = - Y^pi. JdrJ da,[Q{sf{T)~a,)Q{-sf{T)+<j,) 



J={A,B) I •' 

where is given in Eq. ([6|) and the four-vector Lagrange 
multipliers /^'^(t) impose the constraints given in Eq. 

Varying the action ^ with respect to x^ yields the 
usual equation of motion for a string in Minkowski space- 
time (away from the junction), namely the wave equation 



Xi - X, 



where 



K(«,)+^r(«o], (10) 



(11) 



From the conformal gauge conditions ^ the "left" and 
"right" movers satisfy 



a? = 



U2 



0. 



(12) 
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Furthermore, varying the action with respect to Xj, im- 
posing the temporal gauge and using the boundary condi- 
tions gives energy conservation equation at each junction 
(see [4^ and appendix [XI for details): 



0. 



(13) 



As the junction is moving, some of the strings will have 
s/ > while others sf < 0. These represent grow- 
ing/shrinking of the string not only in cr-space, but also 
in real space. The rate of creation of one string must 
balance the disappearance of othcr(s). 

Given an arbitrary initial loop configuration (xi(0, (Tj), 
Xi(0, (Ti)), we aim to solve the above equations in time in 
order to give the full loop evolution and hence {xi{t,ai), 
Xi(t, (7i)). We will now concentrate on junction B. Let 
us first define the "communication route" on the loop for 
a time t at a junction B by 



(14) 



For a time t < junction B receives no influence from 
the dynamics at junction A, with the initial conditions 
specifying the incoming waves at B, so determining 
the behaviour of the junction. For times beyond f^^^, on 
the other hand, vertex B is affected by the dynamics at 
vertex A and the problem becomes more complex. While 
in Appendix IB] we shall obtain analytical results that are 
valid for all times before the two junctions meet, this 
should be considered a special case and in general nu- 
merical methods are required. 

The general formalism is fully described in appendix 
Kl The amplitudes of the outgoing waves are deter- 
mined by the amplitudes of the incoming waves b'^'^ at 
the junction, and we find 



(if + l)ar = &r(l-^f)- 



Efc-"fe 



(15) 



2. Multiple junctions 

Finally, we now generalize the above formalism to a 
string loop with more than two junctions. The impor- 
tant point here is that it is no longer possible to choose 
the (Ji's all to increase (or all to decrease) into a given 
junction. This can be seen, for example, by considering 
carefully the triangular section of the loop shown in figure 
1 . Hence we must treat cases in which there is a different 
orientation between the three strings at a junction. 

To do so we associate a further parameter (5/ with the 3 
strings meeting at junction J. If 5f = +1 then on string 
i, Oi increases into the junction. If, on the other hand, 
6i — —\ for string i then oi decreases into the junction. 
In the action, (5/ therefore ensures invariance under cr- 
reparametrisations, and full details are given in appendix 
El Note also that the values of bf for two junctions that 
are connected by the same piece of string are related. 
The energy conservation equation at the junction is now 
generalized (see appendix [X)) to 



(19) 



For example, consider ^ 62 ~ +1 but changing the 
a orientation of string 3. Then S3 picks up a minus sign 
but (^g = —1, so the energy equation is unchanged and 
energy is still conserved. Finally, the formula for the 
as a function of the incoming waves at a given junction 
J is 



(t) ^sAi 



(Ej-M,)M.(l-cf(t))^ 
M.EfcMfe(l-c((0) 



where 



cf(0 



(20) 



(21) 



with i j ^ k and 



b' 



1 



^ if ^ = -1 



while the evolution of sf is determined by 



(E,M,)M^(1 



'it)) 



'it)) 



where 



cfit) 
Ml 



Kiviit))-K{viit)), 

il^2 - A^3)^ 



l4 



(16) 



(17) 
(18) 



and cyclic permutations. Note that it follows from (|T6|) 
and imposing causality that one must have Mj > 0. This 
is automatically satisfied by (p, q) strings with tensions 
satisfying Eqs. ^ and ([7]). At vertex A the procedure is 
similar, though the incoming waves are now given by the 
af. 



B. Numerical approach 

While the above equations allow for analytical calcula- 
tions for some particularly simple initial conditions (typ- 
ically semi-infinite straight strings) (40l . l4l| , here we con- 
sider more complex loop configurations for which numer- 
ical methods are necessary. We proceed in the following 
way. For every string i connecting two junctions, we work 
entirely with a' and b', reconstructing the closed string 
position x{t,(7i) and velocity x(<, cr,) only when neces- 
sary. The initial conditions fix the initial string tensions 
fii satisfying the triangle inequalities Mj > 0, as well as 
a,^(ai) and h'^{<Ji) between all the junctions. Then we 
calculate the c/ (t = 0) from which {t = 0) is deter- 
mined using Eq. (|^ . Then at time St, sf(5t), uf{6t) 
and u/ {6t) can be calculated. The last step is to use 
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(or the appropriate generalized formula in the case of dif- 
ferent relative orientation of the strings, see for example 
eq. ()A8P ) at each junction to extend the domain of defi- 
nition of a'; (it) and h^(v). The time loop then continues. 

Our simulation ends whenever the length of one string 
goes to zero (so whenever two junctions meet). However, 
the field theory simulations described in the next section 
can of course continue beyond this time. 



C. Initial conditions 

We consider below two different initial conditions that, 
as we shall describe, are closely related. 

The first is the butterfly configuration shown in Fig.[Tl 
consisting of only two junctions (in other words, the small 
triangle shown in the magnifying glass is not present in 
this initial configuration). More specifically, this planar 
configuration consists of a straight string with tension iiq 
and two circular arcs with equal tensions fj,i. We take 
the strings to be initially static at all points: at the site 
of the Y-j unctions this requires Sj = which is satisfied 
when the vector sum of tensions at the junction is zero, 
namely 



j 



(22) 



This corresponds, physically, to the fact that there can 
be no change in momentum occurring within a small vol- 
ume surrounding the junction if the situation is static, 
and is in fact the same condition as derived in Appendix 
A, Eq. (jA3P for the case of the static string. This con- 
dition yields a relationship between the string tensions, 
the radius r of the two arcs making the butterfly "wings" , 
and the distance x of straight string from their centres: 



Mo 
2^1 



X 

r 



R. 



(23) 



The second initial condition is a modification of the 
first initial condition, in that we add, at each junction, 
the triangular shape shown in the magnifying glass in 
Fig. [1] Hence this initial configuration consists of 9 
strings and 6 junctions, all of which are static. All the 
strings within the small triangles are taken to be arcs of 
circles, and are all given the same tension fi2- The rea- 
son for considering this second initial condition is that 
our field theory simulations of the simple butterfly con- 
figuration (with initially only two junctions) will show 
(see section HV C|) that the central straight string can be 
dynamically unstable into splitting into a configuration 
very much like the one shown in figure 1. This dynami- 
cal splitting cannot be accounted for in the Nambu-Goto 
equations, and hence we are forced to put it in by hand 
as an initial perturbation. For that reason the initial size 



of this perturbation is given by a free parameter h, (effec- 
tively the distance between one of the old and newly cre- 
ated junctions), which does not affect the general physi- 
cal behaviour if initially small; this only depends on the 
relative tensions of the strings. The fact that the string 
is initially static sets all angles, shown in Fig. [1] to be 
functions of the tensions iii and the parameter h. 

Finally, we should note that the use of circular seg- 
ments to describe the small perturbation is somewhat 
arbitrary, and it restricts the range of allowed angles (or 
effectively the value of the tensions /i2 in the perturba- 
tion) that can be used to construct this second static ini- 
tial configuration. However, we believe that this choice 
does not have any dramatic effect on the local dynamics, 
since it is only the local curvature at i = around a given 
junction that really affects its evolution. We will illus- 
trate this in section HV C[ when studying the stability of 
Y-j unctions. 



III. FIELD THEORETIC APPROACH 

A. The field theory model 

While the Nambu-Goto formalism is relatively easy 
to analyse numerically, it does not necessarily give a 
complete description of Y-junctions — for instance, one 
might expect important interactions between the strings 
close to and at the junctions, and these are not included 
in the Nambu-Goto action. For that reason we also study 
the butterfiy configuration, listed above, using a field the- 
oretic model of gauge strings that permits the formation 
of junctions. 

Strings with junctions can be formed in numerous dif- 
ferent symmetry breaking schemes: here we study bound 
states and corresponding Y-junctions of the double U(l) 
model [13, which has Lagrangian density: 



(24) 



The notation used follows Refs. [39|, |49[ in that (j) and 
■0 are complex scalar fields, each belonging to one of the 
two Abelian Higgs models that are coupled together via 
the final term. The gauge-covariant derivatives in the 
two halves of the model are: 



and the anti-symmetric field strength tensors are: 



77 

■F iiv 



(25) 
(26) 



(27) 
(28) 
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for gauge fields and i?^ . Finally, 77 and v are constants 
that set the energy-scales of the two halves of the model 
and e.g. \i and k are dimensionless coupling constants. 

For K = the two halves are uncoupled and the topolo- 
gies of the vacuum manifolds, each obeying a local U(l) 
symmetry, require that local string solutions are present 
[l[ (see Refs. for reviews), and take the form of the 

Nielsen-Olesen vertex [52|. For the 0-lialf of the model, 
as a loop enclosing the string is traversed in space, with 
a net non-zero winding of the phase of 4> equal to 27rp 

£ Z), then 101 is forced to zero along the string in or- 
der for the field to be smooth. As a result the field does 
not lie on the vacuum manifold l^l = ?7 near the core and 
the string carries significant potential energy. It addition- 
ally carries a quantized amount of pseudo-magnetic flux, 
equal to 2'Kp/ e in the case of a </> string, since the gauge 
field acts to reduce the phase gradients and so acquires 
a significant curl close to the string. 

For K in the range < k < A1A2, two parallel 
strings from each half can coalesce to reduce the potential 
energy and hence there are bound string states in this 
model [3^. Note, however, that for k > \\/\i\2 the 
potential term is unbounded from below and the model 
hence becomes unphysical. 

B. Collision of straight strings 

The intersection of straight strings has been studied 
in the above model [4§|. It was found that composite 
states do indeed form and that their growth rates at late 
times are well predicted by the Nambu-Goto formalism. 
However, the Nambu-Goto solution differs from the field 
theory case as the moment of collision approaches, (see 
also 150| for a similar analysis in the Type I regime of the 
Abelian Higgs model.) This is partly because any attrac- 
tion of one string to the other just prior to the intersec- 
tion event means that the strings are no longer perfectly 
straight [i^. There are additionally other model depen- 
dent effects involved in the interaction and it is hence 
no surprise that the CKS approach differs from the field 
theory simulation in predicting precisely when the inter- 
section results in composite formation and growth. 



C. Numerical approach 

The numerical approach employed for the field theory 
simulations is largely that of Ref. [i^l , but with a very dif- 
ferent set of initial conditions. That is, we use an exten- 
sion of the Moriarty et al. approach [53| for the Abelian 
Higgs model with the field evolution then preserving the 
analogue of Gauss' law in each half of the model to ma- 
chine precision. We employ a uniformly-spaced cubic lat- 
tice, with spacing Ax, on which the scalar fields are de- 
fined and then the spatial components of the gauge fields 
arc defined on the links between sites, with the time com- 
ponents set to zero via a gauge choice. Energy conscr- 




FIG. 2: The initial (f> phase choice for a planar (1,0) loop, 
ensuring a winding of 2-k in the desired locations. 



vation is accurate so long as the uniform timestep Ai is 
sufficiently small. 

We obtain the initial conditions required for the but- 
terfly conflguration by first setting up the appropriate 
windings in the scalar field and then applying a period of 
dissipative evolution in order to relax the configuration to 
the minimum energy configuration. That is, during this 
period there is an extra term in each of the equations of 
motion that is proportional to the first time derivative of 
the corresponding field and so removes energy from the 
system. Additionally we fix the modulus of the scalar 
field in the region close to the desired centre lines since 
otherwise the configuration would simply contract to a 
point during the dissipative evolution. We apply reflec- 
tive boundary conditions throughout the simulation. 

The initial choice for the phases of the (/) field for a pla- 
nar loop of (1,0) string is made as shown in Fig. [21 If a site 
is above the plane of the loop then it is given the phase 
7r/2, if it is below the plane then it is given 37r/2, while if 
it is in the plane of the loop then it is given either tt if it 
is outside the loop or zero if it is within it. This ensures 
the correct winding structure of the field but it obvi- 
ously yields artificially high gradients on the plane of the 
loop. The modulus of the scalar field is initially chosen so 
that the field lies on the vacuum manifold, except close 
to the string centre lines, as will be explained momen- 
tarily. During the dissipative evolution the gauge field, 
which is set to zero initially, quickly grows to counter 
these phases gradients, while the phase and modulus of 
(f> rapidly adjust themselves in order to minimize the en- 
ergy. Obtaining a (0,1) loop can be achieved by simply 
swapping tp for (p in the above argument, while a (1,1) 
loop is yielded by setting up the phases appropriately in 
both fields. Higher winding numbers cannot be achieved 
by the direct application of the above approach and will 
be discussed below. 

The modulus of the scalar fields is set inside a tube 
around the string centre-line according to the solution 
for an infinite straight string. For a winding 27rm in 
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FIG. 3: The fluxes present in the butterfly configuration 
yielded by the superposition of a (1,1) loop and a (1,-1) loop. 

the phase of (f> and 27rn in the phase of V', this has the 
following form for small displacements r from the string 
centre: 

^{t) « Cr-, (29) 
V'(r) w Dr". (30) 

The constants C and which are dependent on the 
choice of m and n cannot be found analytically, but 
are solved for using essentially the approach of Ref . [s^ . 
Note that if m is finite but n is zero, then even though 
there is no winding in ip, its modulus is still less than 
V near the string as this lowers the total potential term 
energy. However, does remain finite as r — > 0. In 
principle we could fix ji/)! close to the string in this case 
also, but we choose not to since it would not greatly aid 
the fixing of the string position and we wish to minimize 
the artificial restrictions enforced. 

The butterfly configuration illustrated in Fig. [3] can be 
constructed by the superposition of a (1,1) loop and a 
(1,-1) loop after a period of dissipation. Since the equa- 
tions of motion are non-linear there is no precise means 
to do this, however a good approximation is simply to 
sum the gauge fields from each loop and A^^ to give 
the total: 

A, = Al^A-, (31) 

where the + and — refer to each loop. Then for the scalar 
fields: 

^ = ^ ^, (32) 

results in a superposition of complex phases [10, [H, [13, 
[53} . Furthermore, at distances far from any string set 
up in (such that the field is approximately constant 
and close to its vacuum) the form of (p is essentially that 
found in 0_. Using these equations, the time derivatives 
must then superpose as: 

dtA^ = dtAl + dtA;^, (33) 
i^dtP = <P+dtP^ + P-dt(b+. (34) 



K 0^80 0.95 

^(i,o)/27rr7^ 0.864 0.728 

/i(i^i)/27rr;2 1.452 1.133 

m2fi)/2Tiif 1.622 1.271 

7?[(l,0) + (0, 1) (1,1)] 0.840 0.778 

i?[(l, 1) + (1,-1) -» (2,0)] 0.559 0.561 



TABLE I: The energy per unit length of an infinite static 
string, with parameters Ai = A2 = 2e^ = 2g^ — 2, rj = u 



and two values of k as indicated. Note that an Abelian Higgs 
string of unit winding would yield fj. = 2Tvri^ [Hj . R is defined 
in Eq. pSjl. 

While the wings are largely unaffected by this process 
and remain close to the minimum energy solution, a fur- 
ther period of dissipation is required to relax the central 
region, because of the significant interference between the 
two loops. For the case illustrated in Fig. [3] this includes 
the cancellation of the fluxes in t/j along the central string, 
which greatly reduces the energy per unit length of that 
segment. 



IV. RESULTS 

We shall describe our results in units such that the 
initial radius of the circular arcs is unity, since this is a 
useful reference length scale. For the field theory simula- 
tions we set 2 = Ai = A2 = 2e^ = 2q ^, since these values 
arc the standard choices [s^, [13, [49| and neatly yield the 
relevant values of k in the range < k < 1. We have 
studied a number of values of k in this range, with 0.8 
and 0.95 being the two main ones. We additionally set 
rj = V and hence there is complete symmetry between the 
two halves of the model, while we then choose the value 
of rj to set the ratio of the string width (roughly 1/rj) 
relative to the arc radius. Note that although the two 
Abelian Higgs models are each individually in the Bo- 
gomolnyi limit [H^j, with finite coupling between them, 
the -0 field affects the energy per unit length of a (to, 0) 
string even though there is no winding in ip and for ex- 
ample (2,0) strings are stable (see table 

In order to compare field theoretic results with those 
from the Nambu-Goto simulations we additionally cal- 
culate the string tensions in the case of infinite straight 
strings via the method of [1^. We then use these as ^0 
and /.ii in the Nambu-Goto case, although of course we 
are also free to choose arbitrary values for these in the 
Nambu-Goto case. The tensions used are shown in Table 

m 

For moderate values of k (say 0.8) the ratio of com- 
posite to total constituent tensions, i?, for the case of 
(1,0) + (0, 1) (1, 1) is relatively close to unity. Even 
if K is 95% of its maximum physical value then R = 0.78 
and there is only a reduction in energy per unit line of 
about one fifth when a (1,0) string and a (0,1) string 
combine. Further increasing k does not greatly reduce 
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FIG. 4: Results using the Nambu-Goto method with tensions 
set to match a field theory (1,0) + (0, 1) — > (1, 1) case with 
K = 0.8 (left plot), and all tensions equal (right plot). The 
later case corresponds to R = 0.5 and includes a magnified 
region. 



R since most of the energy actually stems from the co- 
variant derivative term. However a large binding energy 
can be achieved if the junction involves flux cancellation 
and therefore we also consider (1,1) + (1,-1) (2,0). 
For K = 0.95 this yields R = 0.56, which indicates that a 
(2,0) is only slightly heavier than a (1,1) string (which 
of course has the same tension as a (1, —1) string). 



A. Direct comparison of field theory and 
Nambu-Goto strings. 

Before presenting our results it is useful to note that 
a circular Nambu-Goto loop of unit radius collapses to a 
point in a time 7r/2. This is also the timescale for collapse 
of the circular arcs in the butterfly case, though strictly 
this is limited to regions of the arcs that remain causally 
disconnected from the junctions during this period (see 
next section). For the case equivalent to (1, 0) + (0, 1) 
(1, 1) with K = 0.8, the Nambu-Goto results shown in 
Fig. m (left) reveal that the central bridge string collapses 
on a timescale shorter than 7r/2. That is if i? = 0.840 
the arcs are still largely intact by the time t = 1.12, 
when the length of the central bridge string is reduced to 
zero. In contrast for the case corresponding to (1,1) + 
(1,-1) — > (2,0) with K = 0.95, the more stable bridge 
remains intact until just prior to the arc collapse: t = 
1.56 rather than 7r/2 « 1.57. A very similar case to this is 
shown in Fig. [3] (right), but R has been set to 0.5 instead 
of 0.56. This is because we cannot strictly continue the 
simulation beyond the time of bridge collapse and this 
choice of R yields a bridge collapse time that is just longer 
than the arc collapse time and hence the approximate 
moment of arc collapse can be seen in the figure. 



In Figs. [5] and [6] we compare the field theory results 
for these two cases with the aforementioned Nambu-Goto 
ones. The snap shots shown are for equally spaced times 
and hence additionally show the velocity of the strings. 
The primary observation is that the agreement between 
the field theoretic and Nambu-Goto results is excellent 
for the times shown, allowing us to extend the results 
from Ref. [49| to strings with global curvature rather than 
just straight strings with kinks, as was previously consid- 
ered. 

It is also instructive to plot the physical length of the 
central bridge string as a function of time, as in Fig. [T] 
For the Nambu-Goto case this can be found analytically 
via the method of Appendix |B] and, since this string is 
stationary, it is just equal to the invariant length: s^—s^. 
Of course, it can also be obtained from our Nambu-Goto 
simulations, and provides a nice test of them, however 
in the field theory simulations, the measurement of the 
string length is non-trivial. 

Fortunately, the geometry of the central bridge is very 
simple since it merely lies along the y-axis at all times. 
Further we can detect the path of the strings by searching 
for lattice grid squares around which there is a net wind- 
ing in the phase of one or both of the scalar fields (j) and ip 
(using the gauge- invariant method of Ref. [IHl)- For the 
case of (1, 0) -I- (0, 1) — *■ (1, 1) we then search for net wind- 
ings that lie within a certain threshold of the y-axis and 
calculate the bridge length as the maximum difference 
in y-coordinatcs. There is a small systematic error de- 
pendent upon the threshold and the angles at which the 
string arcs meet the bridge string, but this is not relevant. 
For the second case involving (1,1) + (1,-1) (2,0), it 
is useful to note that, as shown earlier in Fig. [31 there is 
a (0,1) string that passes around the extremes of the con- 
figuration but not along the bridge. The bridge length 
can hence be taken to be the distance between the two 
points where this string meets the y-axis. 

Results again show excellent agreement between the 
Nambu-Goto and field theoretic approaches, and in fact 
right up to the moment of bridge collapse. However, 
there is a small transient departure initially in the (1, 1)-|- 
(1,-1) (2,0) case, but this is largely due to the ini- 
tial conditions employed and the bridge length measure- 
ment technique. The initial conditions in the field theory 
case force the 4> and string centre lines along very spe- 
cific and idealized paths, but close to the junction these 
are not the natural paths, which are less tightly curved. 
Hence when the strings are released at t = 0, they move 
to follow these more natural paths but overshoot slightly 
before being drawn back. Hence the initial disagreement 
in Fig. [7] is just the (0,1) string attempting to follow a 
smoother and less kinked route across the junction and 
therefore moving outwards from it, but then going too 
far and so undergoing a few low-level oscillations. The 
corresponding oscillations for the other parts of the Y- 
junction, in which the two (1,0) strings that compose the 
(2,0) bridge separate slightly, can be seen in Fig. [HI 
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FIG. 5: The evolution of the butterfly configuration (1,0) + 
(0, 1) (1, 1) with K = 0.8, shown at equally spaced time 
intervals: t = 0.000, 0.267, 0.533, 0.800, 1.067, with larger 
configurations corresponding to earlier times. The field the- 
ory solution is shown as a bitmap, representing the cumu- 
lative projection of its energy density onto the plane, while 
the Nambu-Goto solution is shown as a solid black line. The 
field theory simulation had lattice spacing Aa:: = 0.5/r;, with 
rj = 90. 




FIG. 6: As in Fig. [5] but for n = 0.95 and (1,1) + (1,-1) 
(2,0) . 



B. Collapse of circular arcs in regions causally 
disconnected from the junctions 

As briefly discussed above, an initially static, circular 
Nambu-Goto loop of unit radius, collapses to a point in 
Minkwoski space-time in a time t = 7r/2. In fact, the 
Nambu-Goto solution involves the loop radius undergo- 
ing simple harmonic motion with period 27r. That is, the 
loop collapses to a point but each segment of the loop re- 
emerges on the other side of the loop centre. However, 
the apparent period of the motion is just tt because seg- 
ments from one side of the loop are indistinguishable from 
those on the other side. Not surprisingly, the Nambu- 
Goto action fails to accurately describe a field theoretic 
model once the curvature of the string approaches the 
string width and hence the evolution of the two types of 
string differs just prior to < = 7r/2. 



1 .8| ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 
time 



FIG. 7: The length of the central bridge string as a function of 
time for the analytic Nambu-Goto solution (thin), the numer- 
ical Nambu-Goto results (thick, dashed) and the field theo- 
retic results (crosses). The collection of data with lower bridge 
lengths is for the (1, 0) + (0, 1) (1, 1) case with k — 0.8 while 
higher bridge values correspond to (1,1) -I- (1,-1) — > (2,0) 
with K = 0.95. The initial departure of the field theoretic re- 
sults from the Nambu-Goto ones is due to the measurement 
technique and is a transient effect. 

For the butterfly configuration, by the collapse time, 
information about the presence of the junction will have 
travelled around the arcs an angle 7r/2, which is just equal 
to the invariant length traversed for our chosen initial arc 
radius. Since our initial conditions yield an arc length of 
2{tt — cos~^(i?)) then a length tt — 2cos~^(i?) remains 
unaffected by the junctions, which is between and tt 
for R between and 1. Hence a fraction of the arcs will 
collapse to a point, reaching the speed of light for an 
instant, and yield a sharp kink in the string at a point 
in physical space. This is confirmed in our numerical 
simulations and can be seen in the magnifled region of 
Fig. 21 Note that this is not a cusp in the conventional 
sense, in which the string reaches the speed of light for an 
infinitesimal range of a only. In the butterfly case (with 
unit initial radius) the change in angle is simply equal to 
the invariant length involved, and so would be zero for 
an infinitesimal region. 

In the field theory case, the circular arcs initially fol- 
low the Nambu-Goto dynamics closely, as can be seen in 
Figs. [5] and [SI However, the sharp kink at the moment 
of collapse is not reproduced and rather than the strings 
emerging as an arc of " negative radius" , the energy asso- 
ciated with the arc is dissipated away as radiation. How- 
ever, the string velocity becomes highly relativistic and 
so the string width experiences a signiflcant Lorentz con- 
traction. As this happens, discretisation errors become 
increasingly important, hence this region is even trouble- 
some for the field theory simulations. Since it is not the 
focus of this paper and merely an attribute of our ideal- 





FIG. 8: In the field theory, this diagram shows the possi- 
ble decomposition of a (1, 1) + (1, —1) (2, 0) junction into 
three separate Y-junctions. For the Nambu-Goto simulations 
this can be thought of as the second initial condition (magni- 
fied region of the butterfly configuration shown in Fig. [1} , in 
which the Y-junction has split into three Y-junctions. In this 
case, the angles shown would be determined by the tensions, 
h, and the fact that the junctions are initially static. 



izcd initial conditions, wc shall not dwell on this region 
further. 



C. Stability of Y-junctions to decomposition into 
multiple Y-junctions 

An interesting outcome emerges in the field theoretic 
case of a (1,1) -I- (1,-1) (2,0) junction in that the 
Y-junction itself can decompose as illustrated in Fig. [S] 
Whether or not this occurs will depend upon the tensions 
of the (2,0), (1,1) and (1,0) strings (which cover all ten- 
sions involved due to the symmetry in our parameter 
choice) and upon the orientation of the external strings. 
By coincidence R = fiQ /2/ii is approximately constant for 
(1, 1) + (1, —1) (2, 0) over the range of k that we con- 
sider here (0.8 < k < 0.95). As a result R is effectively 
0.56 over this entire range, essentially fixing both the ini- 
tial orientations and the large-scale dynamics. However, 
the ratio of 7?. = ^i/2^2 (where /i2 is the tension of a (1,0) 
or (0,1) string) varies from 0.86 to 0.78 as k is increased 
over that range. Of course fio/2fi2 must decrease accord- 
ingly since R is almost constant, and hence both types 
of composite string become increasingly stable against 
this decomposition. This can be seen in our results for 
K = 0.8, shown in Fig. [9l and those for k = 0.95 which 
have already been presented in Fig. [SI where in the new 
case the junctions decompose and in the previous case 
they remain stable. 



FIG. 9: Results from a field theoretic simulation for (1, 1) -|- 
(1,-1) — > (2,0) with K = 0.8, showing the decomposition of 
the initial Y-junctions. Snapshots are shown for times t — 0, 
0.36, 0.8 and 1.24; each representing the cummulative pro- 
jection of energy density on the plane. The simulation had 
lattice spacing Ax = 0.5/?7, with rj = 90. 



The final state of the system is no longer two loops for 
the cases when the Y-junctions decompose and instead 
the Y-junctions meet such that the bridge has effectively 
been cut along its length. There are then two separate 
loops in which sit at opposite sides of an elongated tp 
loop to which they are each bound. These inner loops 
then collapse before the outer tp loop. 

Remarkably, we can also analyse this instability of the 
Y-junctions to decomposition in the Nambu-Goto pic- 
ture, by changing our initial conditions to start with the 
decomposed state, but with the three Y-junctions very 
close together. The internal strings between these junc- 
tions are taken to be circular arcs since this approxi- 
mately represents the curving effects of the tensions that 
are applied to the ends of these strings. These are the 
initial conditions shown for time i = in Fig. [TUl although 
the junctions arc so close that they cannot be resolved 
in this plot. However the snapshot for t = 0.50 shows 
that in this case, for which the tensions correspond to 
the K = 0.8 field theory values, the three Y-junctions are 
further separated and approximately circular arcs have 
grown between, which are now readily apparent. That 
is, this and later times indicate that the Nambu-Goto 
evolution also yields the growth of the internal strings 
and so matches the field theoretic results. However, ap- 
plying this procedure to the k = 0.95 case, shows that 
the (1, 1) -I- (1, — 1) — > (2, 0) junction is then stable to col- 
lapse, again agreeing with the field theory results, as is 
clear from Fig. [SI The precise numerical results do differ 
for the unstable case however, since in the field theoretic 
case there is a considerable role of the small-scale micro- 
physics in the breakup of the junction and the breakup 
happens very much faster than is seen in the Nambu- 
Goto case. Additionally, the use of circular arcs here is 
artificial and would not be expected to be wholly accu- 
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rate. However, the approximate matching of the critical 
vahie for k at which the breakup occurs can be consid- 
ered a considerable success for the less computationally 
intensive CKS approach, and crucially, as we show be- 
low, it allows us to make a prediction based purely on 
the Nambu-Goto results as to when a junction will and 
will not be stable to decomposition into many junctions. 

Using the Nambu-Goto approach, the stability of Y- 
junctions can be studied analytically, at least for small 
times. In appendix [Ul we describe this procedure in de- 
tail for junction in Fig. [S] We find that the initial 
perturbation in the decomposed state grows or collapses 
according to the angle p (see Fig. [8]), which purely de- 
pends on the original butterfly tensions, the tension of 
the small arc segments, fi2, and the size of the perturba- 
tion /i, in the following way 



2 



(35) 



where TZ = and, as before, R = Therefore, for 
a given pair of tensions fiQ and fii there is a critical ten- 
sion ^2 = fJ-crit (for a small fixed h), for which p = 0. 
Below and above this critical limit, there are two distinc- 
tive behaviours: one in which the perturbation grows (as 
in Fig. I10[) and one in which it either does not grow sig- 
nificatetly or shrinks, and simply resembles the original 
butterfly (as in the field theory simulation of Fig. [51). To 
illustrate this, we can use the analytic results for Si and 
the evolution of the vertex 'X.a^ as a function of the right 
(and/or left movers), which explicitly is given by (see 
Appendix [C)) 



i-^M,(l-^,)b^. 



(36) 



To map the junction movement in real space, it is useful 
to define the angle 



tan(i^) 



(37) 




FIG. 10: The perturbed butterfly loop corresponding to the 
K — 0.8 field theory case of Fig. [5] using a perturbation param- 
eter h — 0.01 - the instability grows and the loop is unstable. 



P2 big enough, it may even collapse faster than the cen- 
tral bridge does. However, for p < the local curvature 
is such that the strings of the triangular perturbation 
grow in cr-space. In real space, vertex Ai initially moves 
rapidly away from the y-axis and almost along the butter- 
fly wing (string 1), which corresponds to an initial angle 
of (p < IT — tan"^(i?/ (1 — Vl — R^) from the .T-axis (see 
Appendix [C] for details). In figures [TUlandfTTl one can see 
this initial evolution. Later in the evolution (when the 
angle ip reaches tt), the segment A1A2 changes from con- 
vex to concave, and the vertex Ai evolves like any other 
point on the big arc segment, hence moving towards the 
centre of the butterfly wing, as one can see in the last two 
plots of Fig. [ini Therefore, for p < the perturbation 
grows for some time (which depends on how negative p 
initially is), implying the original butterfly Y-junction is 
unstable, leading to the criterion for stability based on 
simply obtaining the value for p. 



which shows the direction with respect to the x-axis for 
t > 0. The explicit treatment of this angle is given in 
appendix O The main result is that there exists a dis- 
continuity in ip when going from p > to p < 0, as 
Fig. [TT] shows [53. This shows how the evolution of the 
splitting of the Y-junction in the original butterfly de- 
pends mainly on the initial local curvature of the strings 
involved. 

When p > (see Fig. [S]), strings 2 and 3 are "com- 
peting" in (T-spacc while the butterfly wing (string 1) is 
not contributing much (note that for small times si = 
to first order in t). After some time and in real space, 
the vertex Ai moves downwards with an initial angle of 
ip > TT + tan~^(\/l — R^ /R) (with the equality in the 
limit of p ^ 0) from the cc-axis, as can be seen in Fig. [TT] 
(see Appendix [C] for details of the calculation) . In this 
case the perturbation does not grow and, for a tension 



V. DISCUSSION 

The role of extended objects in the evolution of the 
Universe has been studied for many years owing to their 
intrinsically interesting properties, their crucial role in 
string/M theory, repeated appearance in GUT models 
and their presence in many condensed matter systems. 
Here we have extended previous studies on the dynam- 
ics of cosmic strings to include the situations envisioned 
in models of cosmic superstrings, namely the existence 
of bound-state strings, and the junctions that join them 
together [Tsj. While analytic studies of the superstrings 
themselves have met with limited success [s^, H^, the 
framework for numerical studies of field theory defects is 
well-developed, and simple models of bound-state strings 
[39! can be used to understand the large-scale properties 
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FIG. 11: Numerical (solid lines) and analytic (dashed red 
lines) evolution of the angle ip, for junction Al, for two cases 
with p > and two with p < 0; all close to the critical 
value p = (blue lines are closer to the critical value). For 
p < (bottom curves) the Y-junction is said to be unstable, 
since vertex A moves along the butterfly wing until it reaches 
180° , and then it starts moving towards the centre of the big 
arc, as shown in Fig. 1101 For p > (top curves), vertex 
A moves downwards, leading to a stable Y-junction. The 
analytic approximations are calculated using R = 0.561 and 
equations HClip and ()C12|) . which are linear truncations (in 
time), and only hold for small times since £^l^{t) ~ h. We 
choose h = 0.01. 

of networks [13, S, S, IZ^I , as well as more detailed anal- 
yses on their individual collisions [i^, [H^ . 

The great utility of Nambu-Goto dynamics is the huge 
reduction in the degrees of freedom compared to field 
theory simulations, allowing for numerical computations 
with far larger dynamic range. However, as is well known, 
the Nambu-Goto dynamics break down as a description 
of gauged cosmic strings when two vortices cross, loops 
contract to a point, or when junctions collide. This is 
well established for the case of the usual Abelian cosmic 
strings, and algorithms have been developed (motivated 
by the field theory results) which nevertheless allow the 
Nambu-Goto equations to be consistently used to model 
the evolution of a network of cosmic strings (see Q for 
details). This same procedure has not yet been estab- 
lished for the case of strings with junctions, and has been 
a key goal of this paper. By comparing the modified 
Nambu-Goto equations for such strings, with the field 
theory evolution for similar strings we have been able 
to establish how well they follow each other, and where 
differences emerge. By concentrating on a simple loop 
configuration, that nevertheless encapsulates important 
dynamics, we have been able to explore various prop- 
erties and directly probe the relationship between these 
two approaches, answering questions such as: What hap- 
pens when junctions collide? What is the effect of loop- 
collapse on the wings of the butterfly loop? What can we 
say about the instability of junctions as the bound states 



try to unzip [i^? 

In terms of the general dynamics, we have seen that 
the evolution of strings with junctions are remarkably 
well modelled by the Nambu-Goto action [i^ [4ll, [sif . 
despite the curvature at the junctions being high. How- 
ever, the field theory approach has opened up a crucial 
new instability that is not present in the usual Nambu 
Goto formalism, namely the break up of a junction into 
three new junctions and the corresponding unzipping of 
the composite string. In other words we still need the 
field theory to understand the outcome of collisions of 
strings and collisions of junctions. 

Examining the field theory of strings with different 
binding energies it became clear that for weakly-bound 
composites the junctions could cause the strings to un- 
zip, care of the new instability. Once this was realized it 
became possible to model this within the Nambu-Goto 
dynamics by introducing the potentially unstable trian- 
gular configuration at the junctions, and solving for its 
evolution. Remarkably, this then allowed us to predict in 
terms of the angle p (or in terms of the string tensions) 
when a junction would be unstable to this decay mode, 
and when it would remain stable, and the results agreed 
well with the field theory simulations. With this done 
there is again excellent agreement between the Nambu- 
Goto simulations and the field theory. The key factor 
here is that given this understanding of the instability, it 
is once again possible to perform, with confidence, large- 
scale cosmological simulations of cosmic supcrstrings us- 
ing these modified Nambu-Goto strings, in the same man- 
ner that it is possible to perform simulations of ordinary 
cosmic strings using Nambu Goto equations. 
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APPENDIX A: EQUATIONS OF MOTION 

Working in the conformal gauge provides an elegant 
way to derive the energy conservation equation at each 
junction. Consider the action given in ^ where we now 
introduce the parameter Si related to the orientation of 
the strings at the junction: 

S = [ dr /"da, [e(sf (r)-aOe(-sf'(r)+a,) 
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EE 



J dr//^.[xf(r,.f(r))-X^(r)] 



where Ai{Bi) is the junction where starts(ends). 
Varying the action (|Aip with respect to Xj we have 



E 



(A2) 



The boundary conditions, on the other hand, are 
Si4J.i {x^ + Sijx^) = f^j (evaluated at each junction) 
which, using Eqs. ([TU]) and (|A2[) gives 



^^-[^■jbecome causally connected. However, since the shortest 
communication route is the central bridge string, and the 
symmetry of the butterfly configuration ensures that the 
only dynamics of the central bridge is for it to change it 
length, then simple analytical results can be obtained for 
later times also, as explained below. 

The initial conditions are a string lying on the y-axis 
(string 0) and two arcs of unit circles (strings 1 and 2) 
in the x — y plane. If is useful to introduce an angle 7 
where, for a static initial configuration, cos 7 = —R. We 
then have: 







A [xT{ui) + six>:{vi)] 



(A3) 



= Y.^'^^^ [(i. + i)ar(uf) + (i-iOe(«/)] 



where 

ui{r) 



(A4) 



Note Eq. (jA3[) reproduces Eq. (QUI for the case of a static 
string Si ~ 0. On imposing the temporal gauge t = t = 
x^(r, (Ti), the /i = component of this equation reduces 
to energy conservation at each junction: 



(A5) 



Now, the constraint that the three strings meet at junc- 
tion B leads to 



2i:^ = [sf + l)aT{uf) + [s, - l)hT{vf). 



(A6) 



Then we proceed depending on the relative orientation 
between the 3 strings in question. To illustrate our 
method, consider the specific example where 5i ~ $2 ~ 1, 
83 = —1 for junction B. That is, the incoming waves at 
junction B are h'^, h'^ and Cg. Eliminating the outgoing 
waves a']^, 02 and 63 between (jA3p and (|A6p gives 

(5^M,)^^=A*3(1 + S3^)4^- E M.e(l-^f)-(A7) 



Eliminating from Eqs. (|A6|) and (|A7p we can solve 
for the unknown outgoing waves. This gives 

- /.(I + sf )a'i = 2^2(1 - - 2a*3(1 + s'Ds!^ 



Xo(t = O,cro) = (0,cro,0), |cro|<sin7, (Bl) 
xi(t = 0, (Ti) = (— cos 7 + cos (Ti, sincTi, 0) , Icil < 7, 

X2(t = 0, 0-2) = (cos 7 - cos (72, sin (T2, 0) , \02\<1- 

If we label the lower vertex as A and the upper one as 
B, then Oi increases towards junction B for all strings. 
Because of symmetry it is sufficient to study one junction, 
say B. Therefore, at i = 



= b'o 
a'l = b'l 
= b'2 



(0,1,0), 

(— sin (Ti, cos (Ti, 0), 

(sin (72, cos (72, 0). 



(B2) 



Since the string is simply stationary and on the y axis 
for all times, then the above t result is valid for all t 
and the ordinarily difficult to handle emitted waves are 
simply the waves set by the initial conditions. When the 
junctions do come into casual contact via string 0, the 
situation is completely unchanged and hence simple ana- 
lytical results can still be obtained. Furthermore, string 
is the shortest communication route, while the results 
below can be used to show that the two longer routes 
are always too long to become relevant before the central 
bridge has collapsed. 

Conservation of energy ([T3]) implies that Rs^ 
and hence after integration we have: 



1 



,s^(i)=sin7--(sf(0-7) 



(B3) 



Thus, it is sufficient to determine the function sf (i). 
In this case we have c\ = C2 = cos(sf — i) and Cq = 



^2^1 — ^ Mi) (1 ~ )t>i (A8) 2cos^(sf — i) — 1 and so it is useful to denote A = t — sf , 
^ J ' whose derivative is A = 1 — s^. Therefore Eq. (flBl ) 



and two similar equations for and bg. Squaring these 
equations and using the gauge conditions (|12p we find 
eq. ((201) with 81=62 = ^ and 6-^ = -1. 



whose derivative is A 
implies; 

A = 



l-i?2 



1 + i? cos A 



Therefore Eq. ([16 



(B4) 



APPENDIX B: ANALYTIC RESULT FOR THE 
BUTTERFLY CONFIGURATION WITH TWO 
VERTICES 



As discussed in section III Al analytical progress can 
be readily made for the period before the two junctions 



This simple separable equation can be integrated to give: 
t sin^ 7 = — cos 7 sin A + A — cos 7 sin 7 -I- 7. (B5) 



Together with the definition of A and Eq. IB3[ we now 
have t, and sf specified as functions of the variable 
A. 
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APPENDIX C: STABILITY OF F-JUNCTIONS 
USING THE NAMBU-GOTO APPROXIMATION 

Using the same small time analysis as in Appendix |b1 
we can study the stability of the initial perturbation in- 
troduced in the butterfly configuration. We will consider 
junction Ai (sec Fig. [8]), however this analysis can be 
equally applied to any other junction. The angles shown 
in Fig. [5] are completely determined by the initial equi- 
librium conditions and hence they are defined in terms 
of the tensions and /i2. The position of junction Ai 
in (T-space is 



sf(0) = ^ = cos-iM^ 



(0) 



2^2 



a, 



(CI) 



where a ~ cos ^ and h is the distance between junc- 
tion A\ (or Ai) and junction A3 in Fig. [H 

We will now analytically demonstrate that the be- 
haviour of the perturbation (i.e. whether it grows or 
collapses) depends on whether the angle p is positive or 
negative. For that we will consider both cases and study 
the behaviour of s^^ and X"^i . We will drop out the 
junction index (Ai) from now on for simplicity. 

Case I: p < 

The initial configuration comprises of three strings 
with tensions /ii , [xi and 

b'i(i = 0,0-1) = (sin (Ti , cos cTi , 0) , 

b'2(i = 0,0-2) = (-sina-2, cos 0-2,0), 

b'3(i = 0,(73) = (-cosa3,sin(73,0). (C2) 

At a later time t the incoming waves at junction A are 

b'i(t, si(t)) = (sin(si(t) - t),cos(si(t) - t),0), 
b'2(i,S2(t)) = (-sin(.S2(i)-t),cos(s2(t)-t),0), 
b'3(t,S3(t)) = (-cos(s3(t)-t),sin(s3(t)-t),0). 

(C3) 

Taylor expanding Si we get Si(t) — Si(0) -I- \it^ -f ... (re- 
member Si = initially), and using the relations between 
the angles we find (to first order in t) 

ci — cos(2Q: + 2t), 
C2 = —cos a, 

C3 = -cos(a + 2t). (C4) 



Now, using equation (|16p . linearising in t and defining 
= cos a ~ -^j^ (which is always less than unity due to 
the triangle inequalities) we find 

1 



-S3 



S2 



27^- 1 



t. (C5) 



Case II: p > 

In this case, the initials conditions can be written as 



b'i(i = 0,0-1) = (siuCTi, cosoi, 0), 
b'2(i = 0,0-2) = (-sino-2,cosCT2,0), 
b'3(i = 0,O3) = (-cos (73, -sin 0-3,0). 

Following the same procedure wc find 



si = 0, 



52 



-S3 



1 + 7^ 



t. 



(C6) 



(C7) 



Having the analytic expressions for s; (which of course 
can be extended further than first order in t) we can easily 
find the analytical expression for X using equations (|A3p 
and (|A6p . to obtain the expression 



(C8) 



In order to study the motion of the vertex in real space, 
we define the angle 



tan(v3) = [ ^ 



(C9) 



Since the expression for this angle is very complicated, 
one can take the limit in which the perturbation size h 
tends to zero, and also consider small deviations from 
the p — case, either with positive or negative p. The 
critical tension p2 which leads to p = is obtained by 
setting ([35|) to zero and solving for /X2, resulting in 



Pcrit 



Ml 



(CIO) 



Therefore, in the limit h 
tion (IC8I) reduces to 



2 cos (cos-i(i?-/i)-7r/2)' 

and p2 ~ Pcrit equa- 



Xy 



R 



-i + 3i?2 + yr~R2 



1 + Vl i?2(l + Vl - i?2 - 2i?Vl - R^) 

(CU) 
for jO < 0, and 



Xy 

Xfr. 



VI - R^ 2 + R^ - 2Vl - 



R 



2R^ 



(C12) 



for p > 0. Notice that, as one should expect, in the 
critical tension limit TZ drops out from the expressions, 
and only R = appears. For both cases (p > and 

p < 0), Xx is initially positive, so it is the y direction 
which changes. For p > 0, the vertex Ai moves with an 
initial angle of 



TT + tan 



R 



(CIS) 
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in the critical limit {fi2 = ficrit), and bigger angles for 
A*2 > fJ-crit', therefore the perturbation does not grow. In 
contrast, for p < 0, the vertex Ai moves away from the 
y-axis, with an initial angle of 

if = 77- tan-1 ( ^ ) , (C14) 



which is practically along the butterfly wing. In this case, 
the junctions separate initially from each other and the 
butterfly configuration is unstable. 
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